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ABSTRACT 

The excess of electrons/positrons observed by the Pamela and ATIC experiments gives rise to a 
noticeable amount of synchrotron and Inverse Compton Scattering (ICS) radiation when the e + e~ 
interact with the Galactic Magnetic Field, and the InterStcllar Radiation Field (ISRF). In particular, 
the ICS signal produced within the WIMP annihilation interpretation of the Pamela/ ATIC excess 
shows already some tension with the EGRET data. On the other hand, 1 yr of Fermi data taking 
will be enough to rule out or confirm this scenario with a high confidence level. The ICS radiation 
produces a peculiar and clean "ICS Haze" feature, as well, which can be used to discriminate between 
the astrophysical and Dark Matter scenarios. This ICS signature is very prominent even several 
degrees away from the galactic center, and it is thus a very robust prediction with respect to the 
choice of the DM profile and the uncertainties in the ISRF. 
PACS: 95.35.+d, 95.85.Bh, 95.85.Pw, 98.70.Vc 

Subject headings: dark matter — gamma rays: observations — cosmic rays — radio continuum: ISM 
— ISM: general — Galaxy: general 



The Pamela and ATIC results have recently raised 
a great interest in the scientific community due to 
the possibility that the observed e + e~ excesses could 
be a signature of the, so-far elusive, particle asso- 
ciated to Dark Matter. The raise in the positron 
fraction above 10 GeV until ^100 GeV seen by Pamela 
(|Adriani et al. 2008al) and the excess of the sum of 
e + and e ~ between ~100 G eV and ^700 GeV seen 
by ATIC (jChang et al. 20081) can be hardly explained 
in a standard Cosmic Ray production scenario and, 
instead, seem to point to a new source of e + and e~. 
Hints of this anomaly were reported also by different 
experiments like HE A T dBarwick et al. 19971) . AMS-01 
(lAguilar et al. 20071: lAlcaraz et al. 2000D and PPB- 
BETS ( |Torii et al. 2008f ). In addition. HESS has re- 
cently presented a measurement of the electron spectrum 
in the range 0.6 < E < 5 TeV (jAharonian et al. 20081 ). 
This anom aly can have a standard astrophysical inter- 
pretation (lAtoian et al. 1995L IZhang and Cheng 2001 



iProfumo 2008llYuksel et al. 20081 iHooper et al. 20091) 
or an exotic one i n volving decayin g Jl/iu et al. 2008 , 
Hisanoet al. 2008bl lYin et al. 20081 IChen et al. 2008 . 
Ibarra and Tranl . lHamaguchi et al. 20081) or the 



annihilation 



Mardon et al 



of 

"20091. 



D M particles (IHisano et al. 2008a 
iZurek 2008L iCholis etaL_2008a 



Bergstrom et al. 2008aL lArkani-Hamed et al. 20091. 
Meade et al. 2009. Ishiwata et al. 2008allHu et al. 20091. 



Nomura and Thaler 20081 



lHall and Hooper 2008 



Barger et al. 2009L IdeBoer 20091 ICholis et al. 2008b 
Fox and Poppitz 20081 ). The latter description, in par 



ticular, seems to favor a DM particle in the TeV range 
and with a thermally averaged annihilation cross section 
(<tav) ~ 10~ 23 cm 3 s _1 . However, this scenario faces sev- 
eral difficulties. A first problem is that, differently from 
the positron ratio, no excess is o bserved by Pamela i n 
the antiproton over proton ratio (jAdriani et al. 2008bl ). 



This means that DM decay/annihilation into hadronic 
channels is mainly forbidde n or at least strongl y sup- 
pressed (jCirelli et al. 2008t iDonato et al. 20081 ), and 
hence one has to resort to models in which only the 
leptonic channels are allowed. The second problem 
is that the annihilation rate required to explain the 
anomaly is about three orders of magnitude above 
the natural expectation of (a\v) ~ 3 x 10~ 26 cm 3 s -1 
for a DM thermal relic which accounts for the cos- 
mological DM abundance. This requires either the 
introduction of large annihilation boost factors from 
the presence of galactic substructure, or some en- 
hancing annihilation mechanism like the Som mcrfcld 
process (jLattanzi and Silk 2008tllbe et al. 200l . 

The fact that hadronic channels have to be suppressed 
to explain the Pamela/ ATIC anomaly implies that only 
few (energetic) photons are produced either if the anni- 
hilation takes place through the fi + fi~ or t + t~ channels 
or in the case of the e + e~ channel through the presence 
of Final State Radiation. With the limited contribution 
of gamma rays accompanying the annihilation process, 
the constraints from gamma observations become thus 
quite weak. Anyway, even though only e + e _ were 
produced in the DM annihilation process, these leptons, 
once in the galactic environment, would interact with 
the Galactic Magnetic Field (GMF) and the Interstellar 
Radiation Field (ISRF). Thus they would lose energy 
producing synchrotron radiation in the radio band and 
Inverse Compton Scattering (ICS) Radiation in the 
gamma band. This secondary radiation thus represents 
a comp lementary observable t o constrain the DM 
signal (iBergstrom et al. 2008E llshiwata et al. 2008bl 
Cholis et al. 20 08c. Nardi et al. 20081 Zh ang et al. 20081 . 
Bertone et al. 20081 iBorriello et al. 2008f ). In the fol- 
lowing we will focus on the synchrotron and ICS 
signals which are expected in the galactic halo. With 
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FlG. 1. — The upper panels show the positron fraction and the total e+e spe ctrum for the CR ba ckground and the DM annihilatio n signal 
compared with the Pamela and A TIC data. A com pilation of previous data (HEA T JBarwick et al. 1997|) and AMS-0 1 JAguilar et al. 20071) for the 
positron fraction and PPB-BETS iTorii et al. 20081) . AMS-01 (lAlcaraz et al. 20001) and HESS dAharonian et al. 20081) for the e + e" spectrum) is also 
shown. The lower right panel reports the gamma spectrum for the CR background and the ICS signal from DM electrons in the Halo together with 
the EGRET measurements and the errors expected after a 1 yr survey by Fermi. The red dashed curve shows the spectrum of gamma-rays produced 
directly through the annihilation into (i~ . The decomposition of the CR background into the IC, bremsstrahlung, pion decay and cxtragalactic 
components is reported as well. The lower left panel shows the DM s ynchrotron emissi on, in units of brightness temperature, 10° away from the 
GC compared with the galactic backgrounds as measured by WMAP IGold et al. 200gf) and the r.m.s. fluctuations o f the CMB. The open point s 
indicate the 9 Planck frequencies while the dotted line shows the expected Planck sensitivity for a 14 months survey (Planck Collaboration 2006). 
The second set of open poi nts indicates the WMAP freq uencies. For comparison it is shown the signal from the WMAP Haze 10° degrees away 
from the GC as derived in ( Doblcr and Finkbcincr 2007). Furthermore wc report the decomposition of the Galactic backgrounds into the dust, 
free-free and synchrotron components together with the synchrotron background derived with Galprop. A model with a WIMP of m x — 1.5 TeV 
which annihilates only into ^ fi~ with a rate (ctav) ~ 5 X 10 
in the text. 

respect to focussing on the Galactic Center (GC) this 
approach provides much more robust predictions due 
to the weaker dependence on the choice of the DM 
profile and thanks to the smaller uncertainties on ISRF 
and GMF. The relevance of ICS signal in relation to 
Pamela has been, i ndeed, discussed in recent papers 
(jZhang et al. 20081 ; iCholis et al. 2008d) which show 
the presence of some tension with the EGRET data as 
well. In the following, we will stress how the situation 
is expected to change with the new data from Fermi 
and, further, we will investigate the peculiar spatial 
distribution which the DM signal is expected to produce. 

We use for the calcula tions a slightly modified ver- 
sion of Galprop v50.1p flStrong and Moskalenko 19981 ; 
iMoskalenko and Strong 1997t ). which solves numerically 
the electron diffusion-loss equation and produces the ICS 
and synchrotron maps. The code also provides maps 
of the CR gamma diffuse emission using available data 
on the CR abundances and the distribution of galactic 
gas. For our calculations wc employ a diffusion coeffi- 
cient D = D (E/E )- a with D a = 5 x 10 28 cnrV 1 , 
-EVj = 3 GeV and a = 0.33, corresponding to a Kol- 
mogorov spectrum of turbulence. The transport equa- 
tion is solved in a cylinder of half-height z = ±4 kpc 



and radius R = 20 kpc, while the GMF used to de- 

1 /2 

rive the synchrotron radiation is modeled as (-B 2 ) = 
Bq exp(— r/rB — \z\ / zb) with Bq = 11/iG re = 10 kpc 
and zb — 2 kpc. It is worth reminding, however, that 
electrons have typically a quite short propagation length 
(in terms of the galactic size) corr esponding to a path 
of 0(1 kpc) (|Delahave et al. 2007f ) before losing a sig- 
nificant percentage of their energy. Thus the final spec- 
trum and distribution of electrons keep only a weak de- 
pendence on the chosen propagation parameters. The 
GMF, on the other hand, is still affected by large uncer- 
tainties esp ecially in the inner kpe's of the galaxy (see 
(|Han 2009f ) and reference therein for a recent review). 
The synchrotron radiation, which is quite dependent on 
the GMF, shares, thus, a similar uncertainty on the nor- 
malization. The InterStellar Radiation Field, which is 
the photon target that determines the ICS signal, is, in- 
stead, better known and the derived ICS signal is thus 
a more robust prediction than the synchrotron signal. 
The IS RF implemented in Galp rop is described in de- 
tails in (jPorter and Strong 2003 ). Finally, for the DM 
profile we choose a very conservative isothermal cored 



one, namely p(r) 
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Fig. 2. — Sky map (in healpix format IIGorski et al. 2003) ) of 
galactic gamma backgrounds at the energy of 10 GeV (top). The 
same with the inclusion of the DM annihilation contribution (cen- 
ter) or Decaying DM (bottom). 

density po = 0.4 GeV/cm 3 at the solar position r Q = 8.5 
kpc. We fix r c = 2.8 kpc for the the core size. However, 
this particular choice is not crucial since we are going to 
calculate the signal not in GC, but the one coming from 
the Halo where the uncertainties on the details of DM 
profile are less relevant. 

We choose to study a single benchmark model with a 
WIMP of m x = 1.5 TeV annihilating in the p + pT chan- 
nel only with a rate (ctav) ~ 5 x 10 _23 cm 3 s~ 1 . The 
resulting electron/positron injection spectrum dN e /dE 
has approximately a constant behavior in energy with a 
cutoff at the mas s of the WIMP. We calculate dN e /dE 
with DarkSUSY (|Gondolo et al. 20041) . which, in turn, 
uses a tabulation of the s pectrum of the decay prod- 
ucts derived with Pythia (jSjostrand 2007| ). The elec- 
tron source term for Galprop is then given by Q(r, E) = 
p 2 (ctav) /2m 2 , x dN e /dE. This model provides a rea- 
sonable good match with the Pamela and ATIC data. 
It is certainly possible to achieve a better fit with a 
mixing of the various leptonic channels, or with particu- 
lar alternative annihilation mechanism or, further, with 
a fine tuning of the propagation parameters 1 . How- 

1 During the review procedure of our paper the Fermi collabo- 
ration has reported a measurement of the e+e — flux in the same 



ever, since the aim of our paper is to focus on the sec- 
ondary radiation, the final results would be only weakly 
affected by the above details on the WIMPs annihila- 
tion process. The results are illustrated in FigfT] The 
upper panels show the comparison of the model with 
the Pamela and ATIC data and with a compilation of 
previous data showing that, indeed, the agreement is 
good. The secondary radiation results are shown in the 
lower panels. The right one shows the expected differ- 
ence between the CR gamma background and the ICS 
produced by the population of DM electrons distributed 
in the galactic halo toget her with the EGRE T mea- 
surements (as taken from ()Strong et al. 2004bf )V Fur- 
thermore the decomposition of the CR background into 
the IC, bremsstrahlung, pion decay and extragalactic 
comp onents is also shown. The extragalactic compo- 
nent (Srccku mar et al. 19981) is from the r eanalysis of 
the EGRET data from ( jStrong et al. 20"04al) . The small 
error bars arc a forecast for T=l yr of data taking 
by Fermi assuming the effective area a s function of en- 
ergy as taken from (|Atwood at al. 2009f ) (roughly A e ff= 
8000 cm 2 above ~1 GeV) a field of view of 2.4 sr and 



no CR contamination hence N~ 



T x fov x /a x 



J AE A eff (E)dN 1 /dE(E)dE. dN 1 /dE(E) is the gamma 
ray flux while /a is the fraction of area of the sky 
where the signal is integrated. The Poisson error is then 



Finally the errors are shown for a logarithmic 



binning of the energy. 

It is worth noticing that the errors expected for one 
year from Fermi survey are tiny enough to detect the ex- 
cess with an high degree of confidence. Even more impor- 
tantly, this excess comes from the halo region, placed sev- 
eral degrees away from the GC and thus in a region were 
the uncertainties on the DM profile are expected to be 
much smaller. Also the u ncertainty on the ISRF, which 
seems anyway not critical (jPorter and Strong 20051 ). nat- 
urally decreases moving away from the GC. A possi- 
ble problem is, in principle, the fact that the DM ex- 
cess can be mistaken with a not well understood CR 
gamma background. Inde ed, the situation is similar to 
the EGRET GeV excess (jHunter et al. 1997j ) which, in 
principle ca n be explained eith er with an "optimized" 
CR mo del (IStrong et al. 2004bl) or with a DM contri- 
bution (jdeBoer et al. 2005I ) 2 . In this case, however, the 
IC excess produced by Pamela/ ATIC is more properly 
a "10-100 GeV excess". Moreover, it generally exceeds 
already the EGRET data, although by an amount which 
is still in principle within the EGRET systematics. A 
more crucial difference is however the spatial distribu- 
tion. While the GeV excess is almost isotropic in the 
sky, the ICS excess has the shape of a circular Haze re- 
flecting the DM distribution in the Halo. This difference, 
indeed, is quite striking, as can be seen clearly in Figj2j 



energy range of ATIC IIAbdo et al. 20091) . The spectrum measured 
by Fermi confirms an excess with respect to the conventional cos- 
mic ray model although the excess is less prominent and smoother 
than the one reported by ATIC. For this broad smooth excess a 
better fit can be achieved through an annihilation into t + t - in- 
stead of /i + fi~ . Using the t+t~ channel as benchmark model, 
however, produces just minor changes in the results derived in the 
following. 

2 Note, anyway, that preliminary results from the Fermi collab- 
oration seem not to confirm the GeV excess. See e.g. the talk 
presented on behalf of the Fermi collaboration at the January 2009 
meeting of the AAS. 
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The CR background instead is expected to lie mostly 
along the galactic plane where the astrophysical sources 
are located. 

The lower left panel shows the DM synchrotron emis- 
sion in units of brightness temperature (T oc v~ 2 F l/ ) 
10° away from the GC compared with the galactic back- 
grounds. We use the WMAP background maps (CMB 
subtracted) and th eir decomposition into synchrotron, 
free- free and dust ([Gold et al. 2008I ) 3 . For illustration 
the frequency spectra in the plot are extrapolated also 
outside the WMAP frequency coverage. We also show for 
comparison the background synchrotron emission calcu- 
lated with Galprop which, indeed, exhibits a close match 
with the WMAP synchrotron spectrum in the 20-100 
GHz range. It has to be noticed that the synchrotron 
galactic CR emission dominates the background only up 
to a frequency of ~ 60 GHz, then there is a small fre- 
quency window which is dominated by free-free (ther- 
mal brcmsstrahlung) emission, while above ~100 GHz 
the background is dominated by dust emission. The 
fluctuations of the CMB dominates around ~100 GHz 
depending on the galactic latitude. The high quality 
data from WMAP, however, allow to efficiently clean 
this further "background". The DM synchrotron radi- 
ation would exhibit in principle a peak with respect to 
the synchrotron b ackground around a frequency ~ 10 5 
GHz( as shown in (|Zhang et al. 20081 )) . where, however, 
the dust background is dominating by many orders of 
magnitude. Restricting the analysis in the more in- 
teresting frequency range < 1000 GHz, the DM sig- 
nal has an almost power law behavior with a slope 
slightly harder than the background, while the spatial 
distribution has a circular shape. These characteris- 
tics in deed correspond to what is found in the WMAP 
Haze dDobler and Finkbeiner 20071 iHooper et al. 20071 ; 
ICumberbatch et al. 20091 ) whose signal we also report 
in the plot for comparison. Notice, however, that the 
Haze feature has still to be firmly established and that 
at the moment it is very much depend ent on the metho d 
employed to separate the foregrounds (jGold et al. 20081 ) . 
Interestingly, we find that, for the GMF model employed, 
the DM signal exceeds the Haze for a factor of ~ 3 simi- 
larly to the IC case. The theoretical signal, on the other 
hand is affected by the uncertainties on the GMF and 
it is difficult to normalize reliably. Moreover, further 
uncertainties come from the systematics involved in the 
separations of the measured signal into the various com- 
ponents, synchrotron, dust, free-free and DM, hence it 
would be difficult to asses the real significance of this 
excess. 

We also consider the case of electrons arising from 
WIMP decay considering a DM signal following linearly 
the halo profile and with the same electron injection spec- 
trum as for the channel. Formally, at the solar po- 
sition, up to diffusion effects, exactly the same positron 
fraction and electron spectrum can be obtained setting 
the DM decay rate to T = po (a\v) /2ni x . The ICS ra- 
diation from the Halo is however significantly reduced 
although Fermi can still discriminate this possibility as 
shown in Fig[2] and Figj3] At this level, however, the 
confusion with a not well understood background could 

3 Data are available at the Lambda web site: 
http: / /lambda. gsfc.nasa.gov/ 
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Fig. 3. — Top panel: Background and DM (either annihilat- 
ing and decaying) latitude gamma profiles averaged in a strip of 
60° along I = compared with the EGRET data. Bottom panel: 
same as above, but with the errors expected with a lyr survey 
from Fermi. At high latitudes the error bars appear artificially to 
increase for the geometry of the 0.5° < |i| < 30.5° strip (which is 
effectively shrinking along b). 



become more problematic although the peculiar circu- 
lar shape of the ICS Haze, present also in this case (see 
Fig©, can help to distinguish the DM signal from the 
astrophysical background. 

Finally, in Fig [3] we report another forecast example 
of the excellent Fermi ability to discriminate among the 
astrophysical and annihilating DM scenario considering 
the latitude profile and a strip of 60° width along 1 = 0. 
We also show in the upper panel the EGRET data in 
the same region and ener gy range (as derived with the 
Galplot package (see also (jStrong et al. 2004bh )). Com- 
pared with the EGRET data the annihilation model 
seems to produce a too much broad peak to fit the data, 
beside producing an excessively high normalization. The 
decaying model is instead difficult to separate from the 
background within the EGRET error bars. With the up- 
coming Fermi data at hands, the analysis can be easily 
generalized to exploit the full angular shape of the IC 
Haze. This would clearly offer the optimal sensitivity to 
disentangle the different scenarios. 

In summary, we have shown that Fermi has the poten- 
tial to test the DM interpretation of Pamela/ ATIC ba- 
sically in a model independent way thanks to the strong 
IC signal which the Pamela/ATIC electrons would them- 
selves produce in the galactic halo. The EGRET data 
seems, indeed, already to disfavor the DM annihilation 
interpretation. Further, the IC signal give rise to a strik- 
ing "IC Haze" feature peaking around 10-100 GeV which 
would provide a further mean to discriminate the DM sig- 
nal from the astrophysical backgrounds and/or to check 
for possible systematics. 
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